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Abstract 

We are interested in some energy functionals concentrated on the discontinuity 
lines of divergence-free 2D vector fields valued in the circle S 1 . This kind of energy 
has been introduced first by P. Aviles and Y. Giga in [3]. They show in particular 
that, with the cubic cost function f(t) = t 3 , this energy is lower semicontinuous. In 
this paper, we construct a counter-example which excludes the lower semicontinuity 
of line energies for cost functions of the form t p with 0 < p < 1. We also show that, 
in this case, the viscosity solution corresponding to a certain convex domain is not 
a minimizer. 


1 Introduction 

1.1 Line energies 

Let 0 be a Lipschitz domain in M 2 . We are interested in measurable vector fields m : fl —» 
M 2 such that 

\m\ = 1 a.e. and V • m = 0 on O, (1) 

where the second equation holds in the distributional sense. In the following, we will 
assume that m is of bounded variation so as to be able to define its jump line. So, we 
consider the set 

A(O) := {m € -BP(0,M 2 ) : \m\ = 1 a.e. and V • m = 0 on 0} . 

Vector fields m E A(fl) are related to solutions of the eikonal equation in 0. Let define 
the set 

S(Sl) := {tp E Lip(Q) : |Vy?| = 1 a.e. and 'S/p E BV( fl)}. 

Then, given m E A(O), there exists a scalar function p E S(Q) such that 

m(x) = ( S/p(x)) ± a.e., 
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where (Vy?) 2- = i?V p stands for the image of T7 p by the rotation R of angle 7r/2 centered 
at the origin in M 2 . Moreover, a function p E Lip(fi) satisfying = m a.e. is unique 

up to a constant and is called stream function. We are now able to define line energies: 


Definition 1 Let f: [0, 2] —>• [0, + 00 ] be a measurable scalar function. Let m E A(Q) C 
.BR(f2,M 2 ). Then, there exists a TL 1 -rectifiable jump line J(m) oriented by a unit normal 
vector u x such that m has traces m±(x) E S 1 on each side of J(m) for TL 1 a.e. x E J{m) 
(see m for more details). 

Then, the energy associated with the jump cost f is denoted by Tf and defined for 
m € A(Ll) as follows: 


Tf(m) = f f(\m + — m-\)d'H 1 (x). 

J J(m) 

f is called the jump cost. Note that the divergence constraint on m E A(Ll) implies that 
for a.e. x E J(m), m±(x) € S 1 and v x satisfy the following condition (see figure [3]): 

m + (x) ■ v x = m-(x) ■ v x . 

Then, in the orthogonal basis (u x , v x ), there exists some angle 6 such that rn± = 
(cos 6, ± sin0) and the jump size is defined as 

t = | m+ — nri-\ = 2| sin 6\. 

Similarly, Tf can be interpreted as a functional of the stream function on the set S(Q): 
Writing m = (Vy?)- 1 € then Tf(m) = Jf(p) where 


V^€5(f2), J f (<p)= [ MV<p)+-(V<p)-\)dH 1 (x). (2) 

An interesting question is to find the minimizing structures of Tf if it exists. Remark 
that for this problem to be relevant, we have to consider a constraint on the boundary 
otherwise all constant functions are minimizers. A natural choice is to minimize Tf along 
all configurations m belonging to the set 


Mo(fl) := {m € M(f2) : m ■ n = 0 a.e. on 50} , 


where n is the exterior unit normal vector of dfl. In terms of the stream function (p, this 
is equivalent to consider the set 


Sq(Q) := {p E S(H) : p = 0 on . 
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1.2 Related models 


The first example of such energy is due to P. Aviles, and Y. Giga. In [5], they have 
conjectured that if /(f) = ^t 3 , then If is the T-limit of the following Ginzburg-Landau 
type energy functional 

. r \ _) [ £|Vu | 2 H—(1 — |rt | 2 ) 2 if u G fP 1 (n,M 2 ) and V • u = 0, 

ALr £ (u) — < Jq e 

I +oo otherwise, 

where D is a bounded open set in R 2 and e > 0 is some parameter. 

For the T-convergence of functionals AG e to If with /(f) = |f 3 , only partial results 
are shown. In [5], the authors have been able to prove the r-liminf property for the 
L 1 convergence using the notion of entropies related to the problem m (see also TO- 
The strong compactness of finite energy sequences has been proved by Ambrosio, De 
Lellis and Mantegazza in [2] and by De Simone, Kohn, Muller and Otto in [7] using a 
compensated compactness method based on a new notion of regular entropy on R 2 . 


The second model we want to address comes from the Ginzburg-Landau theory in thin 
film micromagnetics for some asymptotical regime (see 0 )- Given a bounded domain 
!lcK 2 and a magnetization m = (mi, m 2 , m 3 ): D —>• § 2 , where § 2 stands for the unit 
sphere in R 3 , the energy of m is defined as follows 


E £ (m ) 


e f |Vm | 2 H— 

Jn £ 

+00 



if Vm' = 0 where m' 
otherwise, 


(mi,m 2 ), 


(3) 


where e is a small parameter called exchange length and </ : § 2 —» R is some smooth 
function called anisotropy function. 

If />(m) = |m 3 |“ with 0 < a < 4, only one-dimensional structures are expected and 
it is easy to compute what should be the limiting energy of functionals E e by a 1 D- 
analysis. As for the Modica-Mortola model for phase transition am), E e is expected to 
T-converge to cXj for some c > 0 where /(f) = p = 1 + § is the primitive of y/cj) 
vanishing at 0. The case 1 j>{m) = |m. 3| 2 was studied by R. Ignat and B. Merlet in [9] 
in which a compactness result was proved and a lower bound was found. However, the 
r-liminf property in the definition of T-convergence was established only for limiting ID 
configurations of the form m(x) = ±i' _L for ±x ■ v > 0 with v € S 1 (see figure [3] with 
9 0 = ir/2). 


1.3 Lower semicontinuity, Viscosity solution 

As explained above, some of the line energies If are conjectured to be the T-limit of 
functionals coming from micromagnetics in the space X = L 1 . If that is the case, If has 
to satisfy the following lower semicontinuity property: 
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Definition 2 Let F : X — >• [0,+oo] be a functional defined on some topological space X. 
F is said to be lower semicontinuous or l.s.c. if the following holds: 

V(x n ) n>0 C X, x n — > x => F(x) < liminf F{x n ). 

— n^+oo n—^oo 

Since this property strongly depends on the topology of the space X , we have to 
specify the choice we make for the study of line energies If. 

First of all, due to the non convex constraint \m\ = 1, we need strong compactness 
in L 1 . Moreover, since all the results of the previous part (compactness and r-liminf 
property) holds for the L 1 strong topology, it seems natural to consider the line energies 
If in the space X = L 1 . 

However, since Definition Q] uses the notion of trace of a function, another natural 
choice would be X = BV endowed with the weak topology which is a very common 
choice for phase transition problems. Unfortunately, in the general case, the space BV 
is not adapted to our problem. Suppose /(f) = t p with p > 1 for instance. Then finite 
energy configurations m (i.e. m n — > m in L 1 with If(m n ) < C < Too) are not nec- 

essarily of bounded variation since the total variation of m around its jump line cannot 
be controlled by |m+ — m_| p if p > 1 (see [2]). That is why we need a subspace 

of solutions of the problem CD included in T 1 (H) (and containing BV) because of the 
non convex constraint \m\ = 1 such that we are still able to define a jump line J(m) and 
traces m±. This is done in [6] where a regularity result is shown for solution of CD with 
bounded "entropy production". 

Note that if X and Y are two topological spaces such that Y is continuously embed¬ 
ded in X and F : X —> [0, +oo] is l.s.c. in X then the restriction of F to Y is l.s.c. in Y . 
In this paper, we only want to prove a necessary condition for functionals If to be l.s.c. 
We then prefer to restrict our analysis to HU functions (see remark [2]). 

In the case where /(f) = t p for some p > 0, only partial results are known. In [2], the 
following is conjectured: 

Conjecture 1 Let If be the relaxation of If (only defined on the space BV ) in L 1 : 
If(m) = Inf < liminf If(m n ) : m n € BV and m n —> m in L 1 



If f{t ) = t p with 1 < p < 3 then If is l.s.c. for the strong topology in L 1 . 

For p > 3, this conjecture is false (see ®)- The case p = 3 has been studied by P. 
Aviles and Y. Giga in [5]. More recently the case p = 2 has been proved by R. Ignat 
and B. Merlet in [9]. They also proved that Conjecture [T] holds true for 1 < p < 3 if one 
restricts to configurations m such that the jump size is always lower than \/2. Here we 
are interested in the open case p < 1. 
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We point out that line energies associated with the cost f(t) = t p with 1 < p < 3 
correspond exactly to the expected r-limits of functionals ([3]) when (p(m) = |m 3 |“ with 
0 < a < 4 where Bloch walls seem to be optimal. This is quite natural since when 2D 
structures, as cross tie wall or zigzag wall for instance, have less energy than Bloch walls, 
the T-limit of these functionals may be non lower semicontinuous. In the next part, we 
are going to give a 2D construction which gives some necessary condition on / for If to 
be l.s.c. This condition excludes cost functions of the form f(t) = t v with p < 1; the 
proof is based on a construction in the spirit of [Tj and [12] . 


Theorem 1 Let f: [0, 2] —>• [0, +oo]. Let D be an open and bounded non empty subset of 
R 2 . Assume that If is lower semicontinuous in X = 51/(12, S 1 ) endowed with the weak 
topology. Then f is lower semicontinuous and we have 


f (if) 

limsup - - < 2 limsup f(t). 

t-y 0 t t-y 2 


(4) 


Remark 1 The fact that the lower semicontinuity of If implies the lower semicontinuity 
of f has already been proved in m- The main new point here is the condition Q . 


Remark 2 Theorem [7] is stronger than an equivalent formulation in which BV is re¬ 
placed by some Banach space X such that BV is continuously embedded in X and where 
If is replaced by its relaxation in X. 


As we will see, the lower semicontinuity of functionals If is closely related to the following 
question: Is the viscosity solution a minimizer of If ? In [9j, the authors address the 
following conjecture 

Conjecture 2 Assume that If is l.s.c. in L 1 and that 12 is convex. Then (V^o) -1- is a 
global minimizer of If where ipo(x) = dist(;r, <9I2). 

For a regular domain 12 the distance function <po(x) = dist(x,5J2) belongs to Sq(Q) and 
(V^o) ± is the viscosity solution of the problem ([1]). In particular, if 12 is convex, tpo 
is concave and —D 2 tpo is a positive vectorial radon measure. In [2], the authors give a 
microstructure which shows that the viscosity solution is not a minimizer if f(t) = t p 
with p > 3. As explained below, we are going to give a structure with lower energy than 
the viscosity solution for p < 1. 


Proposition 1 Let f : [0,2] —>• [0, +oo]. There exists a convex domain D such that the 
following holds. Let ipo € 5’o(0) be the distance function (fo(x) = dist(x,3D). Assume 
that (po is a minimizer of Jf defined by (|2|) . Then f satisfies (j4|). 

Corollary 1 There exists a convex domain D such that the viscosity solution is not a 
minimizer of If if f(t) = t p with p € [0,1[. 
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2 Construction of a competitor of the viscosity solution 

In order to obtain the inequality (PT|) . we have to construct a domain 0 on which the jump 
size t = 1771+ — rri- | of the viscosity solution along its singular set is very small. Then, 
we find a competitor whose jump size t is close to the maximal possible value t = 2. In 
other words, we want to substitute small jumps by large ones. 

We will use the polar coordinates (r, 9), r > 0, 9 £ [—7r, 7 t] and we will identify R 2 
and C with the usual bijection. Let D be the unit disk and C be its boundary. 

Let 6q be a fixed angle in ]0, vr/2[ and define the two points A = e l6 ° and A' = e~ lS ° 
on the circle C. Define also Ta (resp. Ta<) the tangent to the circle C at the point A 
(resp. A'). We consider the domain S7 delimited by the large arc {e l6 : |#| > 

Ta and Ta 1 (see figure [Tj). In other words D is the interior of the convex envelope 
of C U {i?} where B = Ta C I Ta'■ Define also Oo = SI fl {|@| < 9q and r > 0} and 

r = on n on 0 = [ab] u [a'b]. 



Figure 1: The domain D and the microstructure m 

We now consider two solutions (po and ip in S'o(D) of the eikonal equation vanishing 
on the boundary: 

• po is the usual distance function: \/x € D, po(x) = dist(x, <9D). 

• p defined by: Mx € D, p(x) = dist(x,<9D U C). 
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A 



Figure 2: Viscosity solution m o on 17 

We also denote by mo = (Vt^o)"* - and m = the corresponding solutions of £[]). 

Then mo,m € Ao(fi). 

We now compute If (mo) and If(m) in order to prove that the function (p has lower 
energy than tpo if f(t) = t p with p < 1. 

Heuristic: The idea is that a small jump along a fixed length is replaced by big jumps 
on a small length : This will reduce the energy for subadditive power costs (i.e. /(f) = t p 
with p < 1) which favor "small jumps". Let us give more details. 

For a small angle 6 q > 0, mo only presents small jumps: mo is C 1 out of segment 
[OB] on which the jump size is |mjj~ — mjj"| =: to = 2sin(f?o)- 

On the contrary, m only presents "big" jumps: i.e. jumps whose size is close to 2. 
The singular set of m consists in 3 different lines : [IB] whose length is equivalent to 0q 
and the two curves C\Cg 0 and 7 e 0 (defined below) on which the jump size tends to 2 and 
the length of these lines is equivalent to 28q. 

As a result, the energy of mo is close to /(2sin#o) while the energy of m is close to 
4#o x /(2). A necessary condition for mo to minimize If is then (see Proposition [T]) 

limsup/(f)/f < 2/(2). 

t-> o 

This excludes subadditive power costs. In the sequel, we are going to make precise 
computations so as to get more informations about the critical angle 8q- 
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Energy of m®: The jump line of mo is the segment [OB] and the traces of mo on each 
side of this line are given by mo,± = — e l 0/ 2 ± e a). j n particular, 


If {mo) = /(2sin Oo)\OB\ = 

COS t>0 

Energy of m: The jump line of m is the union of the 3 curves: 

• C d0 = {e ie : \0\ < 0 O }- 

• 7 e 0 := {z G f2 0 : d(z,C 6o ) = d(z,T)} = {z = re l9 : \0\ < 6 0 , d(z,C) = d(z,d fi)}. 

• The segment [IB] where I = fl [OB\. 

First, let us find a polar equation for the curve 7 g 0 : Given z = re l9 such that |0| < 9q 
and r > 1 we have d(z,Co 0 ) = r — 1 , it remains to compute A := d(z, T). 


Since XI is symmetric with respect to the axis (OB), one can restrict to the case 
M = r e l9 with 0 < 6 < 9o. So A := d(z,T) = \z — P\ where P is the orthogonal 
projection of M = re 10 on the segment [AB] : P should satisfy MP = A OA = Xe l9 ° 
and ~MP ■ aP = 0 . We then compute 

Alf> ■ A.P =mP ■ [IS + (M + mP], 

=^{Xe~ i9 ° {-e i9 ° + re i9 + Xe i9 °)}, 

=A[—1 + r cos(#o — 9) + A]. 

Since MP ■ aP = 0 , this implies A = MP = 1 — rcos(f?o — 9). Then we have z € ^q 0 if 
and only if r — 1 = 1 — r cos(9q — 9) and the polar equation of the curve 'yg 0 is given by 


r{9) 


2 

1 + cos( 0 o - | 0 |) 


-9q<9 < 9 0 


Now, we can compute the energy of m along the curve 7 e 0 : 


• dyl#) = \/r(9) 2 + r'(9) 2 d# where we find r'(9) = -— ^ sin(^o —£) . Introducing 

(1 + cos(#o — 9)Y 

the notation a = 9 q — 9, we obtain 


(1 + cos a)- 


2 \/2(l + cos q) d ^, 

(1 + cos a) z 


4cos(a/2) 

(2 cos 2 (a/ 2)) 2 


So dy reads 


d 7 ( 0 ) = cos 3 (a/2)d9. 



The size of the jump at the point 7 (0) is given by 

i(0) = \m+ -m- 1 = |e* (0o+7r / 2) + e^ d+n/2) \ = |e i(0 °- 0) + 1|. 
Using once again the notation a = 9q — 0, this gives 

t(6) = \J (cos a + l) 2 + sin 2 a = \/ 2(1 + cos a) = 2 cos(a/ 2 ). 

We conclude that the energy of m induced by the jump line 7 e 0 is given by 

r6 ° /[2 cos(a/ 2 )] 


l}(m) = 


'-ft 0 cos 3 (a/ 2 ) 


da. 


The energy concentrated on the arc C 0 O is 

Z}(m) = f{2)'H\C do ) = 26vf(2) 
Finally, we compute the energy on the line [IB]: 

lf(m) = /(2sin 6q)\IB\. 

If the distance function is a minimizer of Zf we should have 

Zf(m) — Zf{rriQ) > 0. 

Now, the preceding equations yields 

Zf(m ) —Zf(mo) =Zj(m) +Zj(m) +Zj(m) -Zf(m 0 ), 
rdo /[2 cos(a/ 2 )] 


/ 


_ 0 O cos 3 (a/ 2 ) 


Since |/B| - |05| = -|OJ| = -r(0) = -- 


da + 20 o /(2) + (\IB\ - \OB\) /(2sin0 o ). 


-, this gives 


cos 2 ($o/2) 

t / \ t , \ f 9 ° /[2cos(a/2)] /(2sin0 O ) 

Z f (m) -Z f (mo) = -3-— —da+ 20 o /(2) - 2 , . ■ 

7 - 0 o cos 3 (a/2) cos 2 (0 o /2) 

Hence, if mo is a minimizer of X/-, the following condition should be satisfied: 


/( 2 sin 0 O ) < 0 o cos 2 ( 0 o / 2 ) 


2 sin 0 i 


0 


< 


sin 0 o 

0o 


1 f e ° /[ 2 cos(a/ 2 )] 


da + /( 2 ) 


0 o Jo cos 3 (a/ 2 ) 

x 2 sup{/(f) : 2cos(0q/2) < t < 2 }. 


sin 0 o cos( 0 o/ 2 ) 

Finally, taking the limsup for 0o —>• 0 in the preceding equation leads to ([4]): 

m 


lim sup ■ 
t^o t 


< 2 lim sup/(t). 
t—>2 


This proves Proposition [T] and corollary |T] follows from the fact that the preceding in¬ 
equality holds false for f(t) = t p with p < 1. Note that in this case, we get something 
more precise than Proposition U) 
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Proposition 2 There exists 8 q g]0, 7t/2 [ only depending on p such that for all 
8 e] — 9q, 9q[, the viscosity solution is not a minimizer of If on Qg where tig is the convex 
set constructed in the previous part (9 being the angle (OB, oA)). 


3 Lower semicontinuity of line energies, proof of Theorem 

CD 

The fact that if If is l.s.c. then / is l.s.c. can be found in [9] (Proposition 1). In this 
section we prove that (U|) is a necessary condition for If to be lower semicontinuous with 
respect to the weak convergence in BV on bounded open subsets of M 2 . 

The key is to use the construction rn E A(Cl) depending on 9q of the first part by 
restriction to Po (See figure EJ). The ID transition defined by ([5]) corresponds to the 
viscosity solution mo of the previous part. Given a small parameter e > 0, it will costs 
less energy to substitute the ID transition around its jump line by the microstructure m 
rescaled at the level e (see figure [4j. 



Figure 3: The vector field m on the left and the lD-transition mo on the right 


We are going to prove Theorem |T] when P = (0,1) x (—1,1). The general case will 
follow easily. Fix 9q g]0,7t/ 2[ and define the ID transition mo for a.e. x\ € (0,1) and 
X 2 € R by 

mo(xi, X 2 ) = m± := (t sin 9q, cos 9q) if ± X 2 > 0. (5) 

Then, let us consider the vector field m = mg 0 of the preceding section restricted to Po 
and define the rescaled and extended vector field m for x± € (0,1) and X 2 E R: 


rh(x i,x 2 ) 


— m((cos$o) 1 ^l,(cos0o) 1 X 2 ) if ((cos 00) 1 3 ; l)(cos0o) 1 ^2)sPo ) 
mo (£ 1 , 3 : 2 ) otherwise. 


10 





Note that rh belongs to j 4(0) and is continuous up to the boundary, rh € C (hlo). Then, 
let n be a positive integer and define m n E yl(Il) by aligning n times the vector field 
rh (see figure n. More precisely, for 0 < i < n and x = (x\,X2) E fl such that 
i/n < x\ < (i + 1 )/n, define 

m n (xi,X2) = fh(nx\ — i,nx 2). 

We have m n (x\ , x'2) = 7710(21,2:2) for \x2\ > 1/n and for all x € 0 , |m n (x)| = 1. Con¬ 
sequently, (m n ) n> 0 converge to mo in L 1 (hl). Moreover, \m n \BV(n) = |m|sy(a) so that 
(m n ) n >o is bounded in BV(VL) and weakly converges to m q. 




n 


Xl 



Figure 4: The microstructure m n 

Since m n is obtained by scaling a fixed structure, it is easy to see that I/(m n ) is 
constant. Indeed, Xf(m n ) = n x 1/n Xf{m) = Xf[m). That is why we obtain the 
following condition: assuming Xf is l.s.c., 

Xf(mo) = /(2 sin#o) < lirninf Xf(m n ) = 

J n—>00 J 

In other words, the viscosity solution costs less energy than the construction mg 0 of the 
preceding part. For this reason, we obtain exactly the same necessary condition (J4]) and 
this ends the proof when hi = ( 0 , 1 ) x (— 1 , 1 ). 

In the general case, let D be an horizontal line such that D n hi 7 ^ 0. Up to a 
translation of hi, we can assume that 0 € D, i.e. D = {x E M 2 : X 2 = 0}. Fix 
9q e] 0 , 7 t/ 2 [ and define uiq E BV(hl,§ 1 ) by: 

for a.e. x E hi, 7710 (^ 15 ^ 2 ) := rn± = (=p sin^o, cos 9$) if ± X 2 > 0 . (6) 

Let g£l and r > 0 and define h l q ^ r = (q, q + r) x (—r, r). n being fixed, we obtain a 
microstructure raj r defined on hlg jT . by rescaling and translating our microstructure m n 
defined on (0,1) x (—1,1). More precisely, for a.e. x E hl 9iT ., 

m q /(xi 1 x 2 ) ■■= m n r Q , -y'j ■ 
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Now, there exists two sequences such that 


OO 


d n n c \J n qi , ri 


i=0 


where for i ^ j, £l quri H fi qjtrj = 0. We can now define the microstructure fh n for a.e. x 
belonging to the whole domain by 



if x € rj for some i € N 
otherwise, 


where rho has been defined in (|6]l . It is clear that rh n - > rho in L 1 (0) and in the 


weak BV sense. Furthemore, 

If (rho) = length (fi n D) If (mo) < lirninf If(rh n ) = length (f2 n D) liminf If(m n ). 

J J n—>-oo 1 n—>-oo 1 

Since Q n D ^ 0 and 12 is bounded, one has +oo > length(0 n D) > 0 which imply the 
thesis. 

Aknowledgment. We want to thank our advisor Radu Ignat for pointing us this 
interesting problem and helpful discussions. 
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